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Abstract 

We show that the scalar leptoquark Yukawa couplings generate a significant lepton flavour vio- 
lation. We compute the light scalar leptoquark contributions to the branching ratios (Br) of the 
lepton flavour violating (LFV) decays £ — > idjij and £ — > £'7 with (i,j = e,/i). We discuss the 
role of the relevant input parameters to these decay rates which are the Yukawa couplings (/w) 
with (a = u,c,t), the light scalar mass Mg 1 and the mixing angle sin2#£Q. We investigate the 
experimental limits from (g — 2)„, /i — e conversion and ir — > ev e , fiu^ to get constraint on the input 
parameter space. We predict that the upper limits on the branching ratios of r — > liljlj can reach 
the experimental current limits. We also show that it is possible to accommodate both r — > liijlj 
and t — > ^7 branching ratios for certain choices of LQ parameters. 

PACS numbers: 13.35.Dx, 13.20.-v,13.35.-r, 14.60.Hi 
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I. INTRODUCTION 



Lepton-flavor violation (LFV), if observed in a future experiment, is an evidence of new 
physics beyond the standard model, because the lepton-flavor number is conserved in the 
standard model. Since the processes are theoretically free from the non perturbative hadronic 
effects they provide accurate predictions for the decay rates and the branching ratios (Br) 
of these processes. Furthermore, they are theoretically rich as they carry considerable infor- 
mation about the free parameters of the used model. On the other hand, the experimental 
work which has been done regarding these decays motivates their theoretical studies. For in- 
stance, experimental prospect for fi — > cy is promising with the recent commencement of the 
MEG experiment which will probe Br(yU — > c-f) ~ 10~ 13 two orders of magnitude beyond the 
current limit. B factories search for the decay mode r — > iiijij at the e + e~ experiment with 
upper limits in the range Br(r — > iiljij) < (2 — 8) x 1CT 8 jl|. Searches for r — > /i/i/i can be 
performed at the Large Hadron Collider (LHC) where r leptons are copiously produced from 
the decays of W, Z, B and D, with anticipated sensitivities to Br(r — > fififl) »s 10~ 8 2j. The 
decay fi — > eee of which there is a strict bound Br(/i — > eee) < 10~ 12 is a strong constraint 
on the parameter space [3]. 

The present experimental upper limits for the branching ratios of £ — > idjij and i — > £'7 
decays are given by 



Br(r iiljij) ~ 



10" 



Br(/i — )■ eee) 
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and 



4m 



(1) 



(2) 



Br(r fj,j) < 6.8 x 10~ 8 , 
Br(r ej) < 1.1 x 10" 7 , 
Br(/i erf) < 1.2 x 10~ n 

Within the SM, the Brs of LFV decays are extremely small. On the other hand, the difference 
between the experimental value of the muon anomalous magnetic moment a M = (g — 2)/2 
and its SM prediction is given by[?-9] 



Aa, 



cxp 



(29.5 ±8.8) x 10 



-10 



(3) 



with a discrepancy of 3.4 a. In spite of the substantial progress in both experimental and 
theoretical sides, the situation is not completely clear yet. However, the possibility that the 



present discrepancy may arise from the errors in the determination of the hadronic leading- 
order contribution to Aa M seems to be unlikely as argued in Ref. flO |. There are many 
attempts, in the literature, to explain this discrepancy through considering new physics 
beyond SM 11- Ijj. 

One of the possibilities for physics beyond the Standard Model is the four-color symmetry 
between quarks and leptons introduced by Pati-Salam [jjj]. The prediction of the existence 
of gauge leptoquarks, which are rather heavy according to the current available data, is a 
direct consequence of this symmetry. 

The current bounds on the leptoquarks production are set by Tevatron, LEP and 
HERA 15[. Tevatron experiments have set limits on the scalar leptoquarks masses Mlq > 



242 GeV. On the other hand, the limits that have been set by LEP and HERA experiments 
are model dependent. The search for these novel particles will be continued at the CERN 
LHC. Preliminary studies at the LHC experiments, ATLAS [if]] and CMS jl7], indicate that 
clear signals can be observed for masses up to 1.2 TeV. 

Our aim in this paper is to analyze the branching ratios for all processes given in Eqs.Q- 
((2]) in the context of the LQ model. These LFV processes are generated at loop level through 
exchanging scalar LQ particles which transmit the lepton flavour mixing from the Yukawa 
couplings to the observed charged lepton sector. Previous studies of such decays were per- 



formed extensively by theorists 



18] . In the present study of these decay channels, the light 



scalar leptoquark effects to £ discussed in detail, namely the contributions of 

the photon and Z boson penguins and box diagrams. Also, we include the predictions for 
£ — > £di£i channels correlated with £ — > £ ,r y rates which are interesting within the framework 
we use. Furthermore, we take into account (g — 2) M , n — e conversion and n — > eu e , fiv^ 
constraints imposed on the input parameter space. This is carried here by considering the 
parametrization introduced in 19| for the case of the £ — > £i£j£j decays. 
The paper is organized as follows: In Section UU we list the relevant terms of the scalar lepto- 
quark Lagrangian to the LFV decays and the analytical expressions of the scalar leptoquark 
contributions to a M and £ — > £'~f decays. The analytical results of the LFV decays £ — > £%£j£j 
will be presented in Sec. III. In Sec. IV, we derive the constraints that can be imposed on some 
leptoquark Yukawa couplings obtained using fi — e conversion. The numerical results for r 
and fi decays will be presented in Sec.V. Finally, Sec. VI will be devoted to the conclusion. 
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II. LEPTOQUARK BASICS 



A. Scalar Leptoquark Interactions 

In this section we list the relevant terms of the scalar leptoquark Lagrangian to our 
LFV decay modes. We consider isosinglet scalar leptoquarks. The effective Lagrangian that 



describes the leptoquark interactions in the mass basis can be written as 



20 



2l|: 



£lq = K (h' ai T kiSR P L + h ai T kjSL P R \ e t S* k +ej(h': j Tl Rik P R + Kf^P^j u c a S k (4) 

- eQ^A^ul - ieQsAnSt&Sk + ieQstmOwZ,,^ S k 

(V) - Q(u-)S 2 w )Pr - Q( u c)S 2 w P L )<, 
V / 



SwCw 

where k — 1,2 are the leptoquark indices, T 3 = —1/2, Q u c = —2/3 are quark's isospin and 
electric charge respectively, Qs = —1/3 is the electric charge of the scalar leptoquarks S k , a 
is up-type quark flavor indices, i, j are lepton flavor indices, cw = cos 6w and $w = sin6V- 
The T k) s L(R) are elements of leptoquark mixing matrix that bring to the mass eigenstate 
basis S k : 

s l = Ts fc Sfc, s R = r kt s R si, (5) 



Here Sl(r) denotes the field associated with the ejPL(R)U^ terms in Clq [20|. Note that in 
the no-mixing case (T = 1), Si^) reduce to Sum which are called chiral leptoquarks as they 
only couple to quarks and leptons in certain chirality structures. Finally, the couplings h 
and h' are 3 by 3 matrices that give rise to various LFV processes and must be subjected 
to the experimental constraints. In this work we do not intend to explore the effects of 
all possible leptoquark interactions. Instead, we try to demonstrate that a simple scalar 
leptoquark model can provide rich and interesting LFV phenomena. 



B. Muon anomalous magnetic moment (g — 2)^ 

The LQ interaction can generate muon anomalous magnetic moment and resolve the 
discrepancy between theoretical and experimental results. The corresponding one-loop dia- 
grams are shown in Fig. []Ja)- QJb) where i = £' = \x. The extra contribution to arising 
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(a) (b) 

FIG. 1: Feynman diagrams contributing to (. — > £'7, S/- denotes the scalar leptoquark with k = 1, 2 
and u c a denotes up-type quark with a = 1, 2, 3. 



from the LQ model due to quark and scalar leptoquark one-loop contribution is given by 



at 2 3 2 

r 1 



a=l fc=l 6 fe 



EE A/f2 



K/JVJ 2 ) (Q( u c)F 2 (x ka ) - QsF x {x ka )) 



m 



Re(tia»Kit r s R ,k r k,s L ) {Q( u c)F 3 {x ka ) - Q s F A {x ka )) 



(6) 



In the above expression, N c = 3, Qs = —1/3, Q u c = —2/3. The kinematic loop functions 
Fi (i = 1, ...,4) depend on the variable x ka = m ( u c)/M| fc , their expressions are given in the 
appendix B. 

Clearly, the use of leptoquark contribution to saturate the deviation shown in Eq. (J3]) leads 
to constraint leptoquark masses Ms k (k=l,2), mixing angle 9lq and the Yukawa couplings 
(h h {,) ) 



c. 1 -»• e 7 

In this subsection, we give the expression for the amplitude of £ — > £'7 which is generated 
by exchange of scalar leptoquark. According to the gauge invariance, the amplitude can be 
written as: 

iM = ieu(p 2 ) (f^rlPl + F^ LR P^\ (i<V9>(Pi)<", (7) 

where e 7 is the polarization vector and q = p\ — p 2 is the momentum transfer. For the 
amplitude of leptoquark exchange at one-loop level, as depicted in Fig. [1] with I ^ £', we 
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have 



F 1 

r 2LR 



3 2 



a mi 



167T 2 M% 
a=l k=l b k 



F 1 

r 2RL 



X {Q(uc)F 2 {x k a) ~ QsFl{x ka )) 

-m iu cj(h ae h'* e ,Tl Rjc T k) s L )(Q(u-)F3(x ka ) - Q s F 4 (x ka )) 
F^ LR (h^h',R^L), 



with x ka = m^ uC ^/Mg k . The branching ratio of £ — > £'7 is given by: 



Br(£ -> £7) 



a.. 



(mf - m 2 ,) 3 



F 7 I 2 



IF 7 I 2 



(8) 
(9) 



(10) 



4r(£) 

In our numerical calculations we analyze the Brs of the decays under consideration by using 
the total decay widths of the decaying leptons T(£). 



in. r -»• 

In this section, we present the analytical results for the LFV r decay into three leptons 
with different flavor within leptoquark model. Next, we give the analytical results relative 
to the branching ratios of t~ — > £^£j£j (the analogous results in the muon sector can be 
obtained by means of a simple generalization.) We perform a complete one-loop calculation 
of the r decay width for all six possible channels, r~ — > fi~fi~fi + , t~ — > e~e~e + , r~ — > 
/i~/i + e~, t~ — )■ e~e + /i~, r~ — > /i~ /i~e + and r" — > ji + e~e~ . The contribution generated by 



the 7-, Z-penguins and box diagrams are presented here separate 



22] 



y. Throughout this section 



we follow closely the notation and thr way of presentation of 

First, we define the amplitude for r~(p) — > £^(pi) £j (p2)£~j (P3) decays as the sum of the 
various contributions, 

«4-(t y £^ £j £j ) ^A-y —penguin A-Z— penguin -^-box- (H) 

In the following subsections, we present the results for these contributions in terms of some 
convenient form factors. 
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<s k 



(c) 



U,, 



(d) 



FIG. 2: Photon (a) and Z-penguin (b) and box (c) Feynman diagrams contributing to t~ — > 
E~£j£j, S k are the scalar leptoquark k = 1,2, u c a are type-up quark with a = 1,2,3. The (d) 
(/U — e) conversion Feynman diagram. 

A. The 7-penguin contributions 

Diagrams in which a photon is exchanged are referred as 7-penguin diagrams and are 
shown in Figs. [2]^a) and [2]^b) when V = 7. The amplitude of T~(p) — > i^(pi) ^7(^2)^/(^3) 
decays can be written as 

iA^-penguin = «(pi) [q 2 l^P L + T*P R ) + im T a ^ {T^ P L + 7?P fl )] u(p) (12) 
3' 
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where q is the photon momentum and e is the electric charge. The photon-penguin amplitude 
has two contributions, one from Fig. [2^a) and the other from Fig. EJ^b) diagrams respectively 



as can be seen from the structure of the form factors, 



rpL,R rp(a)L,R. rp(b)L,R . \ ^ 



(13) 



n 



(a)L _ N c Q( u c) 



3 2 



16vr 2 ^ ^ M q 

a=l k=l °h 



(14) 



T, 



(a)L _ N c Q^ u c) 



3 2 



EE 



16tt 2 ^ ^ M 2 

a=l fc=l ^fc 



rrii 



->(a)i? rp(a)L 



T^(h^ ti,R^ L). 



(15) 
(16) 



and, 



3 2 



T l b)L - --^^^T7T h 'ar h 'ai T \ R ,k T k,S R F & {x ka ) 



16tt 2 ^^M 2 

a=l fe=l b k 



(17) 



(6)i _ N C Q< 



3 2 



EE A/f? 



16vr 2 ^^M 2 

a=l fe=l 5 fe 



^(6)i? _ ^(6)L 



m,, 



+ ^ r ^r^ ifc r fc , Sjs -^F 4 (x fca ) 



T^^T^(h^h',R^L). 



(18) 
(19) 



where Xfc a = m 2 a /M| fe . Note that we have not neglected any of the fermion masses. The 
analytical expressions for the loop functions (i = 1, 6) are given in appendix B. 

B. The Z-penguin contributions 

In addition to the photon penguin diagrams discussed in the previous subsection, there are 
other types of penguin diagrams in which the Z boson is exchanged as shown in Figs. 0(a)- 
[2(b). The amplitude in this case can be written as 



iA 



le 



Z— penguin 2 2 2 

1TI yCu/S 



-u{p 1 ) 1 ^{Z L P L + Z R P R )u{p) 



(20) 



z°w°w 



x u{p 2 )^(g L PL + g R P R )v(p 3 ), 
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As before, the coefficient Z L ^ can be written as a sum of two terms from Feynman diagrams 
in Fig. Eta) and Fig. [2(b): 



Z L ' R = Z^ L,R + Z^ L ' R 



(21) 



where, 



Z (a)L 
z (a)R 
z (b)L 
z (b)R 



3 2 



Yl Yl T72- h 'ar h ai T S R ,k T k,S R 



167T 2 ^ ^ M„ 
a=l k=l b h 

Z {a)L (ti ^h,R^L). 



2C R F 8 (x)-m 2 u C L F 7 (x ka ) 



Nr. 



3 2 



^ M 2 ^'arKi^S R ,k^k,S R 



16yr 2 ^ ^ M 

a=l k=l b k 

Z^ L {ti ^ h,R^ L). 



2Q s tan9 



F 8 (x ka ) 



(22) 

(23) 
(24) 

(25) 



the coefficients Cl(r) and gL(R) denote Z boson coupling to charged leptoquark S and charged 
leptons Il(r), respectively and they are given by 



9l(r) = T Z l{r) - Qem sin 2 6 



Cl(r) — ^3L(fl)(« c ) — Q(u c ) sin 2 9 W , 



(26) 
(27) 



where T 3 l{r) and Q em represent weak isospin and electric charge of Il(r), respectively. The 
loop functions Fj (i=7,8) are presnted in the appendix B. 



C. The box contribution 

The amplitude corresponding to the box-type diagram shown in Fig J2(c) can be expressed 

as, 

iAtox = B^uip^YPMpm^h^PMps)] + B r [u( Pi )YPrHp)]HP2KPrv(P3)] 
+ B^[u( Pl )YPLu(p)}[u(p2h,PRv(p 3 )} + B R [u( Pl )^P R u(p)][u(p 2 )^P L v(p 3 )] 
+ B^[u( Pl )P L u(p)}[u{p 2 )P L u(p)] + B R [u( Pl )P R u(p)}[u{p 2 )P R v(p 3 )] 
+ Biiuip^PLuip^u^a^PLvip^)} 

+ B^uip^PRuip^^a^PRvips)]. (28) 

where 

b l,r = b (c)l,r i = lj _ j4 (29) 
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with, 



(c)L 



(c)L 



B 



(c)L 



B 



(c)L 
4 

B (c)R 



327T 2 

64vr 2 



(30) 



a,a'=l k,k'=l 
3 2 



^ ^ Kr h a'j^k,S R ^k', 



Sl 



hl^h'^rt .,rt ,Dn(ml ,m 2 , ,m% ,m% 



a,a'=l k,k'=l 

N, 



k" 



167T 2 



a,a'=l fc,fc'=l 



x D (m Ua ,m ,m s ,m s ) 



(31) 



(32) 
(33) 
(34) 



Again the loop functions Dq and Do are given in the appendix B. 

By collecting all the formulas, the Branching ratios of r~ — > i^ijij can be written in terms 
of the different form factors as 



Br(r" 



9 ^ 

a mt 



iifr+iTfi a +l[ wi 2 +ii?i a )[8to g 



327rr r 

2{TL T R* + rpLrpR* + h ^ + 



■11 



+ izw+izvr 



2mj 

2(|ZV| 2 + I^V| 2 ) 
1 



t(\B^ + \B^\ 2 ) + U\B L 2 \ 2 + \B^ 



+ ^(l 5 3 L | 2 + l^sT) + \{T?B? + TyBlf + ifBf* + ifBf + h.c) 
- -(T*B^* + T 2 L 5f* + T%B%* + T 2 R i?2 * + h.c) 



(ftf^ + £? + B%Z* L g R + £|^ L + h.c) 



+ - [2{TtZlg L + T^Z* R g R ) + T^Z^ + TfZ^ + h.c] 

+ I [ _ 4(t 2 Xs l + r 2 L ^)-2(r 2 L z^ L + ^ + h.c 



(35) 



where T T is the total decay width of r. All the form factors are real. 
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IV. n-e CONVERSION 



[i — e conversion in the muonic atoms is one of the interesting charged LFV process that 
can occur in many candidates of physics beyond the SM . Accurate calculation of the \x — e 
conversion rate is essential to compare the sensitivity to the LFV interactions in different 
nuclei 23J. In this section, we discuss the constraints that can be imposed on the scalar 



leptoquark couplings using \i — e conversion rate. The dominant contribution to the \i — e 
conversion rate is obtained through considering the tree diagram shown in Fig. EJ^d) which 
leads to the effective Lagrangian 



^eff 



3 2 

EE 

a=l k=l 



Ml 



^h a2 h* ai rl Lk T ktSL {eYPL^){Ual^PLUa) 

- -h a2 h'*j:l R ^ kj s L (ePRii)(u a P R u a ) + (h<r> h',R^ L) 



(36) 



where we have used Fierz transformation for chiral fermions. Pr,l = (1 ±7 5 )/2, u a are light 
and heavy type-up quarks and a matrix is defined by = |[7 M ,7^]- The operators in- 
volving Uaj^Ua, Ua'jsUa, or UaO^Ua do not contribute to the coherent conversion processes 
and thus we can drop them and write 

3 



C 



(u a ) 

eff 



E 

0=1 



(c{# e-fP Rf i + C^l ] e 7 "PL/i) u alil u a 



(c$ eP Lf i + eP R ^ u a u 



(37) 



where we have defined 



a 



(Ua) 
VR 



-h a 2Kl 



k b K 



C SR ~ 7> h ^K\ ^2 TT^^Sr^^Sl 



Mi 



(38) 



Cyl^ and Cgl^ can be obtained by the the exchange h -H- h', R <H> L in Eq. (|38l) . The next 
step for the calculation of \i — e conversion is to match the Lagrangian in E q . (137j) to the 



Lagrangian at the nucleon level. Hence we integrate out the heavy quarks 



24J and so the 
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effective Lagrangian in Eq. (|37|) becomes 



C 



(«) 

eff 



'41 erf Prim + 41 eYPLV) ui. 



u 



(39) 



Then, the effective Lagrangian ( |39l) is matched to the nucleon level Lagrangian 25] 



through the following replacements of the operators 23 



24): 



(40) 



where N represents each nucleon (N = p,n), ipjy are the nucleon fields, and G, f are given 
by 



(u) 



(">P) 



5.1, G 



(u,n) 



4.3 



(41) 



Finally, the Lagrangian at nucleon level can be written as 



where we have introduced the following redefinitions for the vector quantities: 



(42) 



fip) 

U VR 


_ /nf(u) 


f («) 
•/vp 


fi n ) 


/nf(«) 


f («) 


flip) 


_ r<(«) 


f(«) 
JVp 


f( n ) 




Au) 

JVn 



(43) 
(44) 
(45) 
(46) 



while the scalar ones read: 



flip) 

U SR — 


rii u ) 


Q^'P) 


f(n) _ 
U SR - 


ri u ) 


s~i{u,n) 


flip) _ 
°5L _ 




M u >p) 


fin) 
U SL — 


r<i u ) 

°5L 


n iu,n) 



(47) 
(48) 
(49) 
(50) 
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In order to calculate the \x — e conversion amplitude we need to calculate the matrix 
elements of , 4wPn arid ipNl^N of the transition between the initial and the final states of 



nucleus [23l 



24j: 



(A,Z\^ p \A, Z) 


= Zp® 




(A,Z\$ n i; n \A,Z) 


= (A- 


Z)p 


(A,Z\<4> pl %\A,Z) 


= Zp® 




(A,Z\4j nl °4j n \A, Z) 


= (A- 


Z)p 


(A^I^tVatIA z) 


= . 





(n) 



(51) 

where \A,Z) represents the nuclear ground state, with A and Z are the mass and atomic 
number of the isotope respectively, while p^ and p^ are the proton and neutron densities 
respectively. Finally, the p, — e conversion rate is given by 24J: 



m, 



+ 



m. 



(52) 



4 (p^S® + Cfl S&A + AC^ R + 4C$l 
4 (pgls® + Cfl S^A + AC® V® + AC^l V™ 

where V^ N \ S^ N ' are dimensionless integrals representing the overlap of electron and muon 
wave functions weighted by appropriate combinations of protons and neutron densities 23 |. 
For phenomenological applications, it is useful to normalize the conversion rate to the muon 
capture rate through the quantity: 

,(Z,A) 



1 capt 

(Z,A) 



(53) 



The current bounds on B„_ e for Titanium atom and Gold atom obtained by SINDRUM 



collaboration are respectively B^ e {Ti) < 4.3 x 10 



-12 



B^ e (Au) < 7 x 10 



-13 



both 



at 90%CL. The numerical values of V^ N \ and T capt for Titanium and Gold atoms are 
listed In Table [B 

V. NUMERICAL RESULTS AND DISCUSSION 

Let us now proceed to analyse and discuss our numerical results. The quark masses are 
evaluated at the energy scale /x = 300 GeV 28| . which is the typical leptoquark mass scale 
used in this work, 



m t = 161.4 GeV, m c = 0.55 GeV, m u = 11.4 x 10" 3 GeV, 
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(54) 



Nucleus 


S^[m 5 J 2 } 


S^[rnJ 2 ] 


V^[rnJ 2 } 


V^\rnJ 2 } 


^capture [10 $ ] 


Tii 


0.0368 


0.0435 


0.0396 


0.0468 


2.59 


Au^ 7 


0.0614 


0.0918 


0.0974 


0.146 


13.07 



TABLE I: Data taken from Tables I and VIII of [2J. 



while we use the following values for 29 j 

a em = 1/137.0359, M w = 80.45 GeV, M z = 91.1875 GeV. (55) 
We assume as in Ref. 31], that all the couplings h and h! are real and equal to each other 20| . 

h = h' = h*. (56) 
We use leptoquark mass splitting A = 500 GeV in our analysis, where A is defined as 



M| — Mj . Consequently, the remaining parameters in the leptoquark model are the 
mass of the light scalar leptoquark Mg 1 , the mixing angle 6lq, and the couplings h a i (a 
=u, c and t). Also, we assume that the scalar leptoquark may explain the discrepancy Aa M 
between the experimentally measured muon (g — 2) M and its SM prediction Eq. ([3]) and hence 
this condition restricts the possible range of the parameters. We have performed a scan over 
all the input parameters, Ms 1 < 1500 GeV, —1 < sin 9lq < 1 and a em < \h qfJ \ 2 < 1. Then, 
after imposing all the existing constraints arising from n leptonic decays and direct search, 
we select all sets of the input parameters producing the same values for (g — 2)^ at la range 
of data. 

In Fig. [HI we show the allowed regions for different type-up quarks contributions which 
are compatible with ajf* = Aa M at la range of data where the red color (green color) region 
corresponds to top-quark (charm-quark) contribution respectively. As can be seen from the 
left panel of Fig. [31 the dominant contribution is around sin2#£Q ~ 0.7 for both top and 
charm quarks. In addition, we find that, sizeable scalar leptoquark effects to the (g — 2)^ at 
la are obtained for the values of M$ l which satisfy Ms 1 ~ 1 TeV for top quark contribution 
and Ms 1 ~ 400 GeV for charm quark contribution. We note also that it is not possible 
to use the up quark loop contributions alone for LQ, since the couplings |/i U(U | are strongly 
constrained by the tt leptonic decays. It has been found that the LHC has the potential 
to discover light scalar LQ with a mass up to 1.2 TeV and where the Yukawa coupling are 
equal to the electromagnetic coupling j^ . 
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FIG. 3: The allowed regions on the (Ms 1 — \h a)1 \ 2 ) plane (left) and on the (Ms 1 — sin2#£Q) plane 
(right) for top-quark (red color) and charm-quark (green color) contributions, taking into account 
aji = Aa M at la. 

In order to find the constraints on the combination of LQ couplings we require that each 
individual LQ coupling contribution to the branching ratio does not exceed the experimental 
current limits on the Br(£ — > £'7) §2§ and fi — e conversion in nuclei. The latter process is 
used to set the strongest constraints on the product h ull h ue which involve the first generation, 
since, in this case, the process is induced at tree-level. On the other hand, the I — > £'7 decays 
which are induced at one-loop by the photon-penguins, Z-penguins and box diagrams [see 
Figs. E]^a)-|2^c)], allow us to constrain the complementary combinations of the couplings 
involving the second and third generation of quarks, namely h a ih a £/, where a = c, t. The 
fi — e conversion process can be also used to set constraints on the second and third quark 
generations. However, we stress that the bounds from the fi — e conversion are suffering 
from being mo del- dependent due to the non perturbative calculations of the nuclear form 
factors, while the bounds which are obtained from the £ — > £'7 decays are not. Taking into 
account (g — 2) M constraint, and experimental upper limits on the (/i — e)Ti,Au conversion 
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rates [26|, |27j at 90%CL., we obtain the following upper bounds 

h Ufl h ue < 4.38 x 1(T 6 (57) 
for Titanium atom while for Gold atom the bound reads 

h Ufi h ue < 6.25 x 10" 6 (58) 

Clearly, the bounds obtained for both Titanium and Gold atoms are of the same order 
and severely constraint the product of the leptoquark couplings h Ufl h UfJ . Our results are 



consistent with the effective Hamiltonian and approximation used in Ref. (3lJ. The products 
of the couplings h a ih a £i where a=c,t and £, i' = r, /i, e, appear only at the one-loop level 
contribution to the r — > (/i, e) r y and /i — > decays. Therefore, they could be larger if 
they are compared with h u/I h ue ones without violating the experimental upper limits on the 
branching ratios. The bounds obtained on these combinations of couplings in the muon 
sector are given by 

h Cfl h ce < 1.22 x 10" 3 , tights < 5.73 x 1(T 3 (59) 
while for the tau sector they become 

Krh Cil < 4.78 x 1(T 3 , h tT h tfl < 8.13 x 10~ 3 , (60) 
h CT h ce < 7.8 x 10" 1 , h tT h te < 8.0 x 10" 1 , (61) 

In the following studies of the r and /i LFV processes, we will use the parameter space which 
is discussed above. We start by investigating the LFV r and fi decay processes generated 
by the same LQ-scalar interactions as those of (g — 2)^. At one-loop level, the LQ-scalar 
gives contributions to the r~ — > fi~fi~fi + , t~ — > e~e~e + , by means of the so-called 7- and 
Z- penguins and box diagrams. 
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FIG. 4: Scatter plots of (a) Br(r — > fififi), (b) Br(r — > eee) and (c) Br(/i — > eee) as a function of 
light leptoquark mass Ms 1 for top-quark (red color) and charm-quark (green color) contributions. 
The horizontal lines of each plot are the current limits of r and (i LFV decay branching ratios. 
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(c) Br(/U eee) and Br(^ — > e^). for top-quark (red color) and charm-quark (green color) 
contributions The vertical and horizontal lines correspond to the upper limits of r and fi LFV 
decay branching ratios. 
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In Figs HI we present our predictions for the branching ratios of (r — » /i/i/i) (a) and 
(r — y eee) (b) as a function of light LQ mass Mg 1 for top-quark (red color) and charm-quark 
(green color) contributions. These plots have origin in Ms 1 = 300 GeV which roughly corre- 

n 

sponds to the exclusion limit obtained at HERA [15[ for leptoquark masses with couplings 
of electromagnetic strenght. As we can see the main contribution comes from the top-quark 
contribution and can reach 2.47 xl0~ 8 for Br(r — > 3/x) and 5.66 xlO -8 for Br(r — > 3e) which 
are comparable with the present bounds. We find that the main contribution to r — > eee 
and r — > nnjl decays is p roduced from the photon-penguins diagrams which were not taken 



into account in Ref. 



31] . In fact, for large LQ mass (m q <C M^), the photon-penguins are 



proportional to h 2 \og{m q /Ms l )/M'^ i which were known as log enhancement in the litera- 



ture [32J. On the other hand, the naive expectation of Z-penguin and box diagrams leads 
to that they are of orders 0(h 2 m 2 / 'MjJ and O^rUg/M^), respectively. The same con- 
siderations regarding log enhancements hold for the r — > e/i~ fi + and r — > /ie~e + processes. 
However, the uppper limits on the Brs of r — e/i~/i + and r jie~e + could be of C(10 -8 ) 
and the order in size is Br(r — > 3e) > Br(r — > ji~e"e + ) > Br(r — >• e~ ji~ > Br(r — > 3/i). 
Since, r — > e~e~fi + and r — > fi~fi~e + are induced by box diagrams then they are expected 
to be small. On the contrary, since the current bound on the \i — > ej decay imposes very 
strong constraints on the related couplings, the predicted Br(/i — > 3e) is rather too small to 
be observed. 

In Fig. |5l we show the correlations between Br(r — > 3fi) and Br(r — > ^7) in the upper 
left panel, (b) Br(r — > 3e) and Br(r — > ej) in the upper right panel, and (c) Br(/i — > 3e) 
and Br(yU — > ej) in the lower panel. We observe that it is possible to accommodate both 
t — > 3£ and r — > £7 branching ratios for certain choices of LQ parameters. This leads to 
simple correlation like 

Br(j -> ii) Br(fi -> e-y) 

for top-quark contribution, which is in agreement with the ratio expected by the dominance 
of the Penguin- type Figj2]^a)-(b). 
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VI. CONCLUSION 



We have studied the muon anomalous magnetic moment, lepton flavor violating muon 
and tau decays i — > t%tjij and i — > £'7 that are generated by scalar LQ interactions. We 
have found that scalar LQ can explain the discrepancy between the experimental value of 
(g — 2) M and its standard model prediction without any contradictions with the experimental 
bound of LFV tau decay processes. The present experimental limits are used to constrain 
the leptoquark parameter space. We set equal couplings and obtain the upper limits of the 
different product of leptoquark couplings by confronting LFV observable with experimental 
results. Our prediction is that r — > 3/i, 3e, e2/x and r — > /x2e get the leading contributions 
from the so-called photon-penguin diagrams and could be of O(10~ 8 ) which can be accessible 
by the presents experiments and the future linear colliders, such as ILC. On the contrary, 
the current bounds on LFV impose very strong constraints on the Br(/i — > eee) and the 
ratio is too small to be observed in the near future. Hence any observation of LFV processes 
in the charged lepton sector, which are being probed with ever increasing sensitivity, would 
unambiguously point to non-standard interactions. Indeed, such indirect observations taken 
in isolation may not imply much on the exact nature of new physics. But a study of possible 
correlations of its effects on different independently measured charged LFV observable might 
provide a powerful cross-check and lead to identification of new physics through LHC/LFV 
synergy. 
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Appendix A: Constraint form tt — > ev e and ir — > fxu^ decays 



We follow 33|, |34j to constrain leptoquark parameters using pion decay data. Form the 
interactions given in Eq. (J3J), we obtain the effective four- Fermi interaction 



C 



eff 



M 2 

Ok 



(Al) 



M 2 

Ok 



eiP R u a )(d b P R i^] 



By using the Fierz transformation, we can rewrite Eq. llAll) as 

1 



-eff 



2M, 



2 Ki h bj T R,k T kA d L,b7»UL,a) (^jT^Z,,, 



+ 



2M 



s k 

2 h a ih b *F R k r ktL (dL t b u R,a){vL,j€,R,i 



(A2) 



On the other hand, the conventional interaction for the n — > Ivi decay in the SM is given by 



U [^(l-75)/][rfV(l-75H + h.c 



V2 

here \V u d\ is the Cabibbo-Kobayashi-Maskawa (CKM) matrix elements between the con- 
stituent of the pion meson and Gf is the Fermi coupling constant. The ratio Rth of the 
electronic and muonic decay modes is [35] 



R 



tii 



^sm(^ + ~+ eu e ) 



2 \ / 2 2 \ 2 

mi \ / mt — mi x 

7r e 



m fi J V m w ~ m 'fi 

1.2352 ±0.0001) x 10" 



1 + 5 



(A3) 



where 5 is the radiative corrections, Thus the ratio R t h is very sensitive to non standard 
model effects (such as multi-Higges, non-chiral leptoquarks). The experimental value of the 
ratio is [29| 



R, 



exp 



1.2302 ±0.004) x 10" 



(A4) 



The interference between the standard model and LQ model can be expressed by 

1 Re(h ue h': e ) 1 1 Re(/i„X;) 1 



R 



SM- 



-lq — Rth + R 



mi 



th 



m u + m d \^2 G F V ud M 2 m e yfi G F V ud M 2 



^R,k r k,L 



(A5) 
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At 2a level, we get 



rtmin <- m m M 2 j L R,k *>,L < K max (AO ) 



k=1 \ m e ^S k "V lv± S k 



where, 



Rmin = -1-06 x 10~ 8 GeV" 2 , (A7) 
R max = 2.45 x l(T 9 GeV- 2 . (A8) 

The total contribution to Rsai-lq must be smaller than the differences between SM and 
experiment within the allowed error limits. 

Appendix B: One loop functions 

The loop functions used in text are given by 

\2 + 3x — 6x 2 + x 3 + 6x log(x)l 

= 1 12(1-,)* < B1) 

, . [l — 6x + 3x 2 + 2x 3 - 6x 2 log(x)l . , 

F2{X) = - : -; -, (B2) 

v ' 12(1 -x) 4 V ; 

F 3 (x) = 2 ^ x)3 [3 - Ax + x 2 + 2 log(z)] , (B3) 

Fa{x) = 2(1 ^x3 [1 " ^ 2 + 2x log(x)] , (B4) 

FJx) = —j— — -r [16 - 45x + 36x 2 - 7x 3 + 6(2 - 3x) log(x)l (B5) 
36(x — l) 4 

[ — 2 + 9x — 18x 2 + llx 3 - 6x 3 log(x)l 

F « w = 1 m^rr — (B6) 

_ 1 - x + log(x) 

- (x _ 1)2 (B7) 



1 

(x-l) 5 



3 — 4x + x 2 + 4x log(x) — 2x 2 log(x) 



(B8) 
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D (x,y,z,k) 



xlog(x) 



y log(y) 



(x - k)(x - y)(x - z) (y - k){y - x){y - z) 
z\og{z) k\og(k) 



D {y,y,z,k) 



(z — k)(z — x)(z — y) (k — x)(k — y)(k — z) 

' - k\og(k)(y - zf + ( - zk" 2 



D (x,y,k,k) 



D (x,y,z,k) 
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+ (y 2 + z 2 )k- y 2 z) log(y) 
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(k — x) 2 (k — y) 2 {x — y) ^ ^ ^ ^ 
+ (x — y) (k 2 — xy) \og(k) — [k — x)((k — y)(x — y) 
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= x 2 log(x) y 2 log(y) 
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z 2 log(2:) log(A;)A; 2 
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